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ABSTRACT 

In  this  study  reacting  charged  particles  with  intermediate  nucleus  as  a target  ( 67 Zn  , I03Rh,  66Zn  and  I00Ru  ) in 
3j Zn(p  ,n)  f2Ga , w^Rh  ( p , n ) , ^Zn(d ,n)  ® \Ga,  l™Ru(CC,ri)  l^Pd  reactions  as  well  as  proton  with  range  energy 

(4-29. 5MeV)  for  67{Zn(p  ,n)  fnGa  reaction,  proton  with  range  energy  (2. 305-39. 055MeV)  for  l^Rh(p,n)  l^Pd,  deuteron 
with  range  energy  (1.8-15.4MeV)  for  ^ Zn(d,n ) l]Ga  and  alpha  particle  with  range  energy  (10.0-25.0  MeV)  for 
l™Ru  (a,  n)  ™ Pd  are  used  according  to  the  available  experimental  data  of  cross  sections  in  Exfor  library. 

We’ve  calculated  the  cross  section  of  above  mentioned  data  and  results  have  been  obtained  by  using 
(Matlab-8. 34014a)  program. 

The  stopping  powers  have  been  calculated  from  Zeigler  formula  by  using  SRIM-2013  with  the  results  of  cross 
sections  to  calculate  the  neutron-yield  for  reactions  depending  on  67 Zn,  66Zn,  llhRli  and  W0Ru  isotopes  as  targets  for 
reactions,  and  also  comparing  between  cross  section  of  those  reactions  to  choose  the  best  reactions  to  produce  67 Ga  and 
I03Pd  isotope  as  shown  in  figure  (5)  and  (7).  We  also  found  that  comparing  between  neutron  yields  of  mentioned  reactions 
to  choose  the  best  reaction  which  has  high  neutron  yield  as  shown  in  figure  (6)  and  (8). 

KEYWORDS:  Cross  Section  (Excitation  Function),  Stopping  Power,  Neutron  Yield,  Data  Evaluation,  Gallium-67, 
Palladium- 103 

INTRODUCTION 

The  study  of  nuclear  reactions  of  charged  particles  ( p,  d,  a)  with  intermediate  nucleus  as  a target  which  has  a 
high  neutron  yield  due  to  the  advantages  of  the  (p,n),  ( d,n ),  (a,n)  neutron  source  and  to  produce  radioisotopes  as 

Gallium-67  and  Palladium- 103.  They  are  very  important  because  it  provides  a basis  for  a wide  range  of  technical 
applications.  Particularly  using  the  reactions  induced  by  intermediate  energy  protons,  it  is  possible  to  produce  directly 
radionuclides  which  have  been  used  in  medicine  and  industry,  recent  decades,  where  widespread  uses  of  diagnostic  and 
therapeutic  radioisotopes  take  place.  Depending  on  the  type  of  radiation,  the  diagnostic  isotopes  are  classified  into  two 
groups:  [>+,  p emitters  (13N,  isO,  ISF,  62Cu,  6SGa,  etc.)  used  in  Positron  Emission  Tomography  (PET),  and  y-emitters  ( 67 Ga , 
"mTc,  I23I,  etc.)  used  in  Single  Photon  Emission  Computed  Tomography  (SPECT). 

The  total  cross  section  of  such  a production  is  also  important  in  accelerator  technology  from  the  point  of  view  for 
radiation  protection  safety  [1], 
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Gallium-67  is  used  for  tumor  imaging  and  localization  of  inflammatory  lesions  (infections)  [2].  Gallium-67 
behaves  in  the  body  in  a similar  way  to  ferric  iron.  It  is  commonly  used  as  a trivalent  citrate  compound  for  nuclear 
medicine  imaging,  and  is  a valuable  agent  in  the  detection  and  localization  of  certain  neoplasms  and  inflammatory  lesions 
[3],  Gallium-67  with  a half  life  of  (T i/2=3.217  d),  its  decays  to  stable  67Zn  by  electron  capture.  Its  decay  emissions  include 
gamma  rays  of  93.3  keV  (37.0%),  184.6  keV  (20.4%)  and  300.2  keV  (16.6%).  is  a cyclotron-produced  radioisotope  for 
which  considerable  demand  exists  [3].  This  radioisotope  67Ga  is  well  known  and  widely  used  in  the  field  of  nuclear 
medicine.  67Ga  has  become  one  of  the  most  frequently  employed  cyclotron  produced  radioisotope  over  the  last  two  decades 
and  is  a widely  used  single  photon  marker  for  detecting  the  presence  of  malignancy  and  for  diagnosing  inflammatory 
lesions  [4]. 

Palladium- 103  is  used  to  make  sealed  seeds  (brae hy therapy)  as  a permanent  implant  seeds  for  early  stage  prostate 
cancer  [2].  Palladium-103  is  extensively  used  in  the  treatment  of  prostate  cancer  as  well  as  in  ocular  melanoma,  and  is 
mostly  applied  in  the  form  of  sealed  seeds  (brachytherapy)  [4].  Palladium-103  has  a suitable  half-life  of  (Ti/2=16.991  d),  its 
decays  to  stable  103Rh  by  electron  capture  100  %.  which  de-excites  by  means  of  a heavily  converted  internal  transition,  as  a 
result  of  both  processes  (EC  and  IT),  X-rays  and  Auger  electrons  are  both  emitted  which  are  ideally  suited  for 
brachytherapy  (its  electron  capture  decay  resulting  in  abundant  emission  of  Auger  electrons  and  low  energy  X-rays  (20-22 
keV))  [3,5], 

NUCLEAR  REACTION 

The  (7-value  of  the  X(a,b)Y  reaction  is  defined  as  the  difference  between  the  final  and  initial  kinetic  energies  and 
is  given  by  [6,7,8]: 

Q = [M„+Mx-(Mi+Mr)]x  931.5  (1) 


The  threshold  energy  for  a nuclear  reaction  is  defined  as  the  smallest  value  of  a projectile’s  kinetic  energy  at 
which  a nuclear  reaction  can  take  place  [8,  9]. 
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NUCLEAR  CROSS  SECTIONS 
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To  study  nuclear  reactions,  it  is  necessary  to  have  a quantitative  measure  of  the  probability  of  a given  nuclear 
reaction.  In  which  the  cross  section  of  nuclides  is  the  effective  area  presented  [10].  When  an  accelerated  charged  particle 
interacts  with  a target  nucleus  a nuclear  reaction  takes  place,  ultimately  leading  to  a stable  or  radioactive  product  nucleus. 
A nuclear  reaction  is  characterized  by  a cross-section,  describing  the  probability  that  a particle  interacts  [9,  11].  The  cross 
section  has  the  units  of  area  as  it’s  the  square  of  nuclear  radius.  The  standard  unit  for  measuring  a nuclear  cross  section  (o) 
is  the  barn  (b),  which  is  equal  to  (lbarn=10-24  cm2)  or  ( lbarn=l(T28  m2  ).  The  reaction  cross  section  data  provides 
information  of  fundamental  importance  in  the  study  of  nuclear  systems.  The  cross  section  is  defined  by  [12]: 
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Were  R is  the  number  of  reactions  per  unit  time  per  nucleus.  I is  the  number  of  incident  particles  per  unit  time  per 
unit  area.  In  general,  a given  bombarding  particle  and  target  can  react  in  a variety  of  ways  producing  a variety  of  light 
reaction  products  per  unit  time.  The  total  cross  section  is  than  defined  as  [8]: 


a 


total 


(4) 


Where  n , is  the  partial  cross  section  for  the  process. 

STOPPING  POWER 


The  energy  of  charged  particle  loses  per  unit  path  length  of  the  material  it  traverses.  Generally  speaking,  any 
charged  particle  can  have  either  electronic,  nuclear,  or  gravitational  interaction  with  the  particles  of  the  material  through 
which  it  passes.  However  the  gravitational  interaction  is  too  low  to  be  of  any  significance  and  is  generally  ignored.  The 
total  stopping  power  is  then  just  the  sum  of  the  stopping  powers  due  to  electronic  and  nuclear  interactions  [13,14,15]. 
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Where  Seiectro„iC  is  the  electronic  stopping  power,  S„uciear  is  the  nuclear  stopping  power.  The  electronic  stopping 
power  is  always  much  larger  than  the  nuclear  stopping  power.  For  most  practical  purposes  the  nuclear  component  of  the 
stopping  power  can  also  be  ignored  as  it  is  generally  only  a fraction  of  the  total  stopping  power.  For  particles  such  as 
electrons,  this  statement  is  always  valid  since  they  are  not  affected  at  all  by  the  strong  nuclear  force. 

Electronic  Stopping 

Electronic  Stopping  Power  (Se)  The  beam  particles  hitting  a target  get  slowed  down  by  interactions  with  the 
electrons  until  they  are  in  thermal  equilibrium  with  their  surroundings.  As  a consequence  one  gets  a wide  neutron  spectrum 
if  the  projectiles  are  completely  stopped  in  the  target,  even  if  all  nuclear  reactions  were  two-body  reactions.  The  loss  of 
kinetic  energy  in  a nuclear  encounter  would  be  much  larger  [7].  In  the  scope  of  this  work,  the  electronic  stopping  powers 
were  calculated  using  the  Ziegler  formulae  [15].  Reliable  data  are  available  for  many  elements  over  a wide  range  of 
energies.  However,  in  order  to  obtain  values  for  all  elements  over  a continuous  range  of  proton  energies,  the  authors  fit 
curves  through  the  available  experimental  data  to  generate  coefficients  for  use  in  a semi-empirical  parameterization  of  the 
stopping  power  as  a function  of  proton  energy  E (keV)  and  the  target  atomic  number  Z2.  Se  is  assumed  to  be  proportional  to 
E045  for  E < 25  keV,  except  for  Z2<  6,  where  it  is  proportional  to  E° 25,  for  25  keV  < E < 10  MeV  [16,17], 

= + (6) 
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S -A  FAl 

° Low  ~ 
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(7) 


Sffigh  (High  energy  stopping)  is 

A5  \n(—  + A1E) 
v E 

High  \.  (8) 


Where  E (KeV)  energy  of  proton  and  the  coefficients  A,  for  each  Z2,  available  from  SRIM  [16,  18] 
For  10  MeV  < E < 2 GeV, 


„ . . ME.  ME^2  a / E . 

S High  — Ag  + A10 ( e ) + Aj ! ( e ) + A12 (9) 

For  use  in  the  parameterization  Proton  stopping,  St('H),  in  various  target  materials  for  10  kev  < E < 100  MeV, 
High-energy  proton  stopping  Se  ~St,  Proton  stopping  power  for  1 kev  < E < 10  GeV  and  Z2<  92,  based  on  a combination 
of  theoretical  calculations  and  experimental  data  (for  20  keV  < E < 1 MeV)  [18,19]. 

Nuclear  Stopping  Power  (Sn) 

A beam  of  charged  particles  bombarding  the  neutral  atoms  of  a target  interacts  with  the  atomic  nuclei  and  atomic 
electrons  of  the  target.  The  ratio  of  the  energy  lost  in  interaction  with  the  atomic  electrons,  to  the  energy  lost  in  the 
interaction  with  the  atomic  nuclei  [14] 

2mp/mes  4xl03 

Thus  the  energy  lost  by  interaction  with  the  nuclei  is  negligible  compared  with  that  lost  by  interaction  with  the 
electrons,  for  solid  targets  with  a given  thickness  there  will  be  the  more  nuclear  interactions  and  hence  the  electronic 
stopping  cross  sections  is  comparable  with  the  nuclear  cross  section  [15].  the  theory  of  ion-solid  interactions  is  well 
established  and  it  has  been  discussed  in  detail  [18,20].  Where  the  authors  also  described  the  computer  program  SRIM 
[16,18].  which  can  be  used  for  calculating  stopping  power  and  range  for  ion/target  combinations  with  Z<92. 

For  sufficiently  high  projectile  energies,  S,  ~ Se.  Sn,  in  units  of  eV/1015  atoms/cm2,  for  any  projectile  with  energy  E 
(kev),  is  given  by  [16,18]. 


S„(E) 


8.462  ZtZ2M  1Sn(e) 
(AT,  + MjKZ,023  + Z2°-23) 


(10) 


Where  the  reduced  ion  energy,  s,  is  defined  as 


Reduced  Ion  Energy  - £ - 


32.53 M,M2(£/M1) 
Z1Z2(Ml  +M2)(Z1°'23  +Z2°'23) 


(11) 


Where  Mi  and  M2  are  the  projectile  and  target  masses  (amu),  and  Zi  and  Z2  are  the  projectile  and  target  atomic 
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numbers.  For  e < 30  keV, 


Sn(£) 


ln(l  + 1.1383g) 

2{e  + 0.013216*0'21226  + 0. 195936*05) 


For  e > 30  keV,  unscreened  nuclear  stopping  is  used,  and  S„  (e)  simplifies  to 


S.(e) 


\n£ 

l£ 


(12) 


(13) 


NEUTRON  YIELDS 

The  neutron  yield  (Yn)  detected  per  incident  particle  (proton,  deuteron  alpha),  for  an  ideal,  thin,  and  uniform 
target  and  mono-energetic  particles  beam  of  incident  energy  Eb  is  given  by  [13]. 

Yn={Ndx)a{Eb)n{Eb)  <„> 


Where  Nd  is  the  a real  number  density  of  target  atoms,  x is  the  target  thickness,  c is  the  reaction  cross  section  and 
r|  is  neutron  detection  efficiency.  If  the  target  is  sufficiently  thick.  For  a target  which  is  not  infinitesimally  thin,  the  beam 
loses  energy  as  it  passes  through  the  target,  if  the  target  is  sufficiently  thick,  and  there  exist  one  atom  per  each  molecule 
(i.e.,/=  1)  and  taking  (E')  = 1,  then  the  resulting  yield  is  called  the  thick-target  yield  which  is  given  by  [13,21,22]. 


Y„(Et) 


° {E  > dE 

clE  / dx 


(15) 


Where  Ethr  is  the  reaction  threshold  energy,  o(E)  is  the  cross  section,  dE/dX  is  the  incident  particle  initial  energy, 
N is  the  atomic  number  of  target  per  unit  volume,  which  is  defined  as  follows: 


N 


w pN 
A 


(16) 


Where,  w is  the  abundant  in  the  combination,  is  the  combination  density,  A is  the  mass  number,  NA  is  the 
Avogadro’s  number. 

Thus,  by  measuring  the  yield  at  two  closely  spaced  energies  (El)  and  (E2),  one  can  determine  the  average  value 
of  the  integrand  over  this  energy  interval  as  follows  [23]: 


y{e2)-y(ex)  = 


°(E) 

dE  / dx 


(£2  - El) 


J e„ 


(17) 


Where  (Eb)  is  the  average  of  (El)  and  (E2).  If  a (E)  are  available  in  the  literature  as  a function  of  projectile  energy 
(Eb)  for  natural  elements,  then  the  neutron  yield  can  be  calculated  using  eq.(13).  If  neutron  yield  is  available  as  a function 
of  projectile  energy  (Eb),  then  eq.  (13)  can  be  used  to  calculate  a(E)  as  a function  of  (Eb).  Thus,  consequently  one  can 
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calculated  the  neutron  yield  by  using  eq.  (17),  For  natural  elements  and  if  only  one  stable  isotope  is  available  in  nature, 
then  [24] 

= Y (E  ) m 

Where  (Yo)  is  the  neutron  yield  per  106  bombarding  particle  for  the  natural  element. 

RESULTS  AND  DISCUSSIONS 

The  3i^72(/7,7z)  32G'<7  Reaction 


The  cross  sections  of  the  Zn(p,n ) f2Gci  reaction  have  been  published  as  a function  of  proton  energy  by 

Szelecsenyi  [25],  Hermanne  [26],  Levkovskij  [27],  Little  [28]  in  Exfor  library  Accordingly  we’ve  used  the  range  of 
energies  above  the  threshold  energy  (1.8105MeV)  with  range  of  energy  (4-29. 5MeV)  in  step  of  (0.5  MeV),  the  cross 
section  of  '',\Zn  ( p , n ) f2Ga  reaction  depends  on  the  kinetic  energy  of  the  proton  and  directly  proportional  with  the  kinetic 

energy  of  proton  until  reaches  to  (11.0  MeV)  which  is  the  maximum  value  of  the  cross  section  (c>Max=739.09  mb)  after  that 
it  decreases  with  the  increasing  of  kinetic  energy  of  proton.  The  aim  of  study  is  for  the  cross  section  and  neutron  yield  for 
reaction  which  is  shown  in  figure  (1)  and  (6)  because  this  reaction  is  one  of  the  important  interactions  that  lead  to  the 
production  of  radioactive  isotope  such  as  67Ga. 


Figure  1:  The  Cross  Sections  of  3 \Zn(p , n)  ^Ga 


The  30  2/1  (d  , ll)  3 [ Gd  Reaction 


The  cross  sections  of  the  ^ Zn  (cl , n ) ” Ga  reaction  have  been  published  as  a function  of  deuteron  energy  by 
Williams  [29]  Accordingly  the  cross  section  is  used  within  range  of  energy  (1.8-15.4MeV)  in  step  of  (0.4  MeV),  the  cross 
section  of  ^ Zn  (d  ,n)  ” Ga  reaction  depends  on  the  kinetic  energy  of  the  deuteron  and  directly  proportional  with  the 

kinetic  energy  of  deuteron  until  reaches  to  (8.2MeV)  which  is  the  maximum  value  of  the  cross  section  (aMax=422.22mb) 
after  that  decreases  with  the  increasing  of  kinetic  energy  of  deuteron.  The  aim  of  study  is  for  the  cross  section  and  neutron 
yield  of  reaction  which  is  shown  in  figure  (2)  and  (6)  because  this  reaction  is  one  of  the  important  interactions  that  lead  to 
the  production  of  nuclei  radioactive  such  as  67Ga. 
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Figure  2:  The  Cross  Sections  of  3*Z«  ( d , n ) ftGa 
The  ( p , n)  U^P d Reaction 

The  cross  sections  of  the  l^Rh  {p  , ft)  Pd  reaction  have  been  published  as  a function  of  proton  energy  by 

Sudar  [30],  Hermanne  [31],  Harper  [32]  , Johnson  [33]  in  Exfor  library.  Accordingly  we’ve  used  the  range  of  energies 
above  the  threshold  energy  (1.3385MeV)  with  range  of  energy  (2. 305-39. 055MeV)  in  step  of  (0.75MeV),  the  cross  section 
of  (p , n)  '^Pd  reaction  depends  on  the  kinetic  energy  of  the  proton  and  directly  proportional  with  the  kinetic 

energy  of  proton  until  reaches  to  (9.805MeV)  which  is  the  maximum  value  of  the  cross  section  (oMax=509.93  mb)  after  that 
it  decreases  with  the  increasing  of  kinetic  energy  of  proton.  The  aim  of  study  is  for  the  cross  section  and  neutron  yield  of 
reaction  which  is  shown  in  figure  (3)  and  (8)  because  this  reaction  is  one  of  the  important  interactions  that  lead  to  the 
production  of  radioactive  isotope  such  as  Palladium- 103. 
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Figure  3:  The  Cross  Sections  of  45  Rh  ( p , Tl)  d 


The  4ARu  (a,n)  46Pd  Reaction 

The  cross  sections  of  the  'j'j'/fft  (d ,n)  'fi'Pd  reaction  have  been  published  as  a function  of  proton  energy  by 
Qaim  [34],  Hussain  [35]  in  Exfor  library.  Accordingly  we  ’ve  used  the  range  of  energies  above  the  threshold  energy 
(7.6815MeV)  with  range  of  energy  (1.8-15.4MeV)  in  step  of  (0.4  MeV),the  cross  section  of  l44Ru(CC ,n)  ^ Pd  reaction 

depends  on  the  kinetic  energy  of  the  alpha  particle  and  directly  proportional  with  the  kinetic  energy  of  proton  until  reaches 
to  (19MeV)  which  is  the  maximum  value  of  the  cross  section  (aMax=526.55mb)  after  that  it  decreases  with  the  increasing  of 
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kinetic  energy  of  alpha  particle.  The  aim  of  study  is  for  the  cross  section  and  neutron  yield  of  reaction  which  is  shown  in 
figure  (4)  and  (8)  because  this  reaction  is  one  of  the  important  interactions  that  lead  to  the  production  of  radioactive  isotope 
such  as  Palladium- 103. 


CONCLUSIONS 

From  the  results  of  the  cross  sections  and  neutron  yield  for  (p,n),  (d,n),  (alpha,n)  reactions  with  many  intermediate 
nucleus  target  and  production  of  many  radioisotopes  with  short  half  life  we  have  been  observed  many  conclusions  for  this 
study  of  (p.n),  (d,n),  (alpha,n)  reactions. 

• The  main  contribution  process  to  produce  67 Ga  radionuclide  [ (7\Zn  ( p , n)  %Ga  and  “ Zn(d,n ) ^ Ga  reactions] 

have  a high  cross  section  value,  the  maximum  value  of  the  cross  section  (om.ix=739.09  mb)  at  kinetic  energy  of 
proton  (11.0  MeV)  for  (p.n)  reaction,  while  the  maximum  value  of  the  cross  section  (oMax= 422.22  mb)  at  kinetic 
energy  of  deuteron  (8.2MeV)  for  (d,n)  reaction  illustrates  in  figure  (5).  the  (p,n)  reaction  is  the  best  to  produce 
Gallium-67  because  of  the  high  cross  section  value  compare  with  (d,n)  reaction. 
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Figure  5:  The  Cross  Sections  For  ^ Zn  (p  ,n)  ^ Ga  and  “Z/i  (d , n ) ffia  Reactions 
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In  figure  (6)  the  ^Zn  ( p , n ) ^Gtl  reaction  has  neutron  yield  value  a better  than  “ Zn(d , n ) ^Ga  reaction 

because  of  the  high  cross  section,  lower  stopping  power  of  incident  proton  on  nucleus  target  67 Zn  than  incident  deuteron  on 
nucleus  target  66Zn  with  dependence  on  equation  (17)  as  shown  in  figure  (5),  The  (p,  n)  as  neutron  source  with 
intermediate  mass  nuclei  as  a target  is  the  best  which  has  the  high  neutron  yield  compare  with  (d,n)  reaction. 
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Figure  6:  The  Neutron  Yield  for  3j Zfl  (p  ,Yl)  32  Gf/  and  j0Zn(d  ,Yl)  3 1 CjCI  Reactions 

The  main  contribution  process  to  production  103 Pd  radionuclide  \'°3Rh(p,n)l^Pd,  '^Ru(a ,n)1^Pd  reactions] 

have  a high  cross  section  value.  The  maximum  value  of  the  cross  section  (c>Max=509.93  mb)  at  kinetic  energy  of  proton 
(9.805MeV)  for  (p,n)  reaction,  and  the  maximum  value  of  the  cross  section  (GMax=526.55mb)  at  kinetic  energy  of  alpha 
particle  (19MeV)  for  (a,n)  reaction  illustrates  In  figure  (7).  the  (p,n),(a,n)  reactions  are  best  reactions  to  produce 
Palladium-103  but  the  (p,n)  reaction  is  better  than  (a,n)  reaction  to  produce  palladium-103  because  it  has  a high  value  of 
cross  section  at  low  energy  (9.805MeV)  while  the  (a,n)  reaction  has  a high  value  of  cross  section  as  well  at  energy 
(19.9MeV). 


Figure  7:  The  Cross  Sections  for  ^Rh(p ,n)  and  l™Ru(a ,n)  l^Pd  Reactions 
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For  '°*Rh  (p  ,n)  '°63 Pd  i ™Ru  (a  ,n)  Pd  reactions,  the  neutron  yield  of  (p,n)  reaction  is  higher  value  than  (a,n) 

reaction  as  neutron  source  as  it  has  a high  cross  section  and  low  stopping  power  of  incident  proton,  the  neutron  yield  value 
dependence  on  equation  (17)  as  illustrates  in  figure  (8),  The  (p,  n)  as  neutron  source  with  intermediate  mass  nuclei  as  a 
target  is  the  best  which  has  the  high  neutron  yield  compare  with  (a,n)  reaction.. 


Figure  8:  The  Neutron  Yield  for  l^Rh(p  ,n)  Pd  and  '^Ru  (CC,n)  l^Pd  Reactions 
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